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We present a measurement of the time-dependent CP- violating (CPV) asymmetries in — > Kg-K" 
decays based on 124 million T(45') — » BB decays collected with the BABAR detector at the PEP-II 
asymmetric-energy B Factory at SLAC. In a sample containing 122 ± 16 signal decays, we obtain 
the magnitude of the direct CPV asymmetry C^o^o = 0.401q ± 0.09 and the magnitude of the 

CPV asymmetry in the interference between mixing and decay Sj^o^o = 0.481q|7 ± 0.06 where the 
first error is statistical and the second systematic. 

PACS numbers: 13.25.Hw, 13.25.-k, 14.40.Nd 



The BABAR 0| and Belle collaborations recently 
reported observation of CP violation in B meson de- 
cays through measurements of the time-dependent CP- 
violating (CPV) asymmetry in B'^ decays into charmo- 
nium final states. In the framework of the Standard 
Model (SM), where CP violation is a consequence of the 
presence of a complex phase in the Cabibbo-Kobayashi- 
Maskawa (CKM) quark mixing matrix these mea- 
surements determine the parameter sin 2/3, with /3 = 
a,rg{—VcdV*^/VtdV^l). The consistency of the observed 
value of sin 2/3 with the Standard Model expectations 
provides strong evidence that the CKM mechanism is 
the dominant source of CP violation in the quark sector. 
A major goal of the experimental studies of B decays 
is to provide additional information to examine the va- 
lidity of this conclusion and search for evidence of new 
physics (NP) in possible deviations from the SM. One 
avenue for the observation of NP is provided by CP vio- 
lation studies of decays dominated by penguin loop-level 
b sqq {q = {d, s}) transitions |3j Sill- While in 
the SM the time-dependent CPV asymmetries in these 
decays measure sin 2/3, additional radiative loop contri- 
butions from NP processes may alter this expectation. 
Presently, the B factory experiments have explored time- 
dependent CPV asymmetries in three such decays, which 
in the SM are dominated by the penguin b sss tran- 
sition: B° r/K° HQ, B° K+K-K° 0, and 
B° (j)K° H The latter results hint at a possible 
deviation from the SM, but are inconclusive. 

In this letter we present the first measurement of the 
time-dependent CPV asymmetries in the decay B'^ 
KgTT^, which has a measured branching fraction B{B'-* 
K°TT°) = (11.9 ± 1.5) • 10-6 [13. The CKM and color 
suppression of the tree-level b — *■ suu transition leads to 
the expectation that this decay is dominated by a top 
quark mediated b —^ sdd penguin diagram, which carries 
a weak phase arg(VtbV'(*). If other contributions, such 
as the b — s- suu tree amplitude, are ignored, the time- 
dependent CPV asymmetry is governed by sin 2/3. The 
deviation from sin 2/3 due to standard model contribu- 
tions with a different weak phase is estimated to be at 
most 0.2 [13. 

The results presented here are based on 124 mil- 
lion T{AS) BB decays collected in 1999-2003 with 
the BABAR detector at the PEP-II e+e" collider, lo- 
cated at the Stanford Linear Accelerator Center. The 
BABAR detector, which is fully described in 0, pro- 



vides charged particle tracking through a combination 
of a five-layer double-sided silicon micro-strip detector 
(SVT) and a 40-layer central drift chamber (DCH), both 
operating in a 1.5 T magnetic field in order to pro- 
vide momentum measurements. Charged kaon and pion 
identification is achieved through measurements of par- 
ticle energy- loss (dE/dx) in the tracking system and 
Cherenkov cone angle {9c) in a detector of internally 
reflected Cherenkov light (DIRC). A segmented CsI(Tl) 
electromagnetic calorimeter (EMC) provides photon de- 
tection and electron identification. Finally, the instru- 
mented fiux return (IFR) of the magnet allows discrimi- 
nation of muons from pions. 

We search for B° -f^s'"'" decays in hadronic events, 
which are selected based on charged particle multiplicity 
and event topology 13]. We reconstruct tt+tt^ 
candidates from pairs of oppositely charged tracks. The 
two-track combinations must form a vertex with 7r"'"7r- 
invariant mass within 3.5cr of the nominal mass [l^ 
and reconstructed proper lifetime greater than five times 
its uncertainty. We form tt" — > 77 candidates from 
pairs of photon candidates in the EMC that are iso- 
lated from any charged tracks, carry a minimum en- 
ergy of 30 MeV, and possess the expected lateral shower 
shapes. Finally, we construct B° — > K'^tt^ candidates 
by combining and tt'^ candidates in the event. For 
each B candidate two nearly independent kinematic vari- 
ables are computed, namely the energy-substituted mass 
TOES = y/{s/2 + PiPB)^/Ef +p%, and the energy dif- 
ference AE = Eg — Vs/2. Here, {Ei,pi) is the four- 
vector of the initial e+e~ system, ^/s = Ef — pf is the 
center-of-mass energy, ps is the reconstructed momen- 
tum of the B° candidate and Eg is its energy calculated 
in the e'^e~ rest frame. For signal decays, the toes dis- 
tribution peaks near the B'^ mass with a resolution of 
~ 3.1 MeV/c^ and the AE distribution peaks near zero 
with a resolution of ~ 40 MeV. Both the toes and the AE 
distribution exhibit a low-side tail from energy leakage 
out of the EMC. We select candidates within the window 
5.2 < TOES < 5.29 GeV/c^ and -150 < AE < 150 MeV, 
which includes the signal peak and a "sideband" region 
for background characterization. For the 1.7% of events 
with more than one candidate we select the combination 
with the smallest — J2i=-jT0 k° ("^» ~ ''^i)'^ /'^rm' where 
rrii {m[) is the measured (nominal) mass and arm the 
estimated uncertainty on the mass of particle i. 

For each B'' ^s7r° candidate we examine the re- 
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maining tracks and neutral candidates in the event to 
determine if the other B meson, -Btag, decayed as a 
or a B'^ (flavor tag). Time-dependent CPV asymmetries 
are determined by reconstructing the distribution of the 
difference of the proper decay times, At = tcp — ttag, 
where the tcp refers to the signal and itag to the 
other B. At the T{AS) resonance, the At distribution 
follows 



' BO 



(At) 



\At\/T 



4t 



(1) 



[1 ± ( S"/ sin {AnidAt) ~ Cf cos (AnidAt) ) ] , 

where the upper (lower) sign corresponds to i?tag decay- 
ing as -B° {B'^), t is the -B° lifetime averaged over the 
two mass eigenstates, Ato^ is the mixing frequency, Cf 
is the magnitude of direct CPV in the decay to final state 
/ and S the magnitude of CPV in the interference be- 
tween mixing and decay. For the case of pure penguin 
dominance, we expect Sx^nO = sin 2/?, and C^o^o = 0. 

We extract the CPV parameters from an unbinned 
maximum-likelihood fit to kinematic, event shape, fla- 
vor tag, and time structure variables. We verified that 
the selected observables are sufficiently independent that 
we can construct the likelihood from the product of one 
dimensional probability density functions (PDFs). The 
PDFs for signal events are parameterized from either 
more copious fully-reconstructed B decays in data or 
from simulated samples. For background PDFs we select 
the functional form from data in the sideband regions of 
the other observables where backgrounds dominate. We 
include these regions in the fitted sample and simulta- 
neously extract the parameters of the background PDFs 
along with the CPV measurements 

The sample of B^ 
by random K'^tt'^ combinations from e+e^ ~* <N i<l = 
{u, d, s, c}) fragmentation. Monte Carlo studies show 
that contributions from other B meson decays can be 
neglected. We exploit topological observables to discrim- 
inate the jet-like e"'"e~ — > qq events from the more uni- 
formly distributed BB events. In the T{AS) rest frame 
we compute the angle 9g between the sphericity axis [isj 
of the B^ candidate and that of the remaining parti- 
cles in the event. While |cos6'g| is highly peaked near 
1 for e^e~ — > qq events, it is nearly uniformly dis- 
tributed for BB. We require jcosfljl < 0.8, eliminat- 
ing 83% of the background. In addition, we include in 
the fit a Fisher discriminant variable, which is defined as 
JT = O.53-O.6OL0+I.27L2, where Lj = ^. |p*|| cos6'*p', 
p* is the momentum of particle i and 9* is the angle be- 
tween p* and the sphericity axis of the B^ candidate. 

We use a neural network (NN) to determine the fla- 
vor of the i?tag meson from kinematic and particle iden- 
tification information 16]. Each event is assigned to 
one of five mutually exclusive tagging categories, de- 
signed to combine flavor tags with similar performance 
and At resolution. We parameterize the performance of 



-ftTgTT^ candidates is dominated 



this algorithm in a data sample (i?flav) of fully recon- 
structed B'^ £''^*^^7r+//3+/aj^ decays. The average 
effective tagging efficiency obtained from this sample is 
Q = J2c 4(1 - 2^"'^)^ = 0.288 ± 0.005, where and w= 
are the efficiencies and mistag probabilities, respectively, 
for events tagged in category c. For the background the 
fraction of events (e^) and the asymmetry in the rate of 
B^ versus B^ tags in each tagging category are extracted 
from the fit to the data. 

We compute the proper time difference At from the 
known boost of the e"*"e^ system and the measured 
Az = zcp — zti^g, the difference of the reconstructed decay 
vertex positions of the B^ ^s""*^ and Stag candidate 
along the boost direction (z). A description of the in- 
clusive reconstruction of the -Btag vertex is given in 'l3|. 
For the B'^ ^s""" decay, where no charged particles 
are present at the decay vertex, we exploit the fact that 
the flight distance of the B meson transverse to the beam 
direction (^ 30 /im) is small compared to the flight length 
along the beam (^ 260 /im). We then determine the de- 
cay point from the intersection of the K'^ trajectory with 
the interaction region by constraining the B vertex to the 
interaction point (IP) in the transverse plane. The posi- 
tion and size of the interaction region are determined on 
a run-by-run basis from the spatial distribution of ver- 
tices from two-track events. The uncertainty in the IP 
position, which follows from the size of the interaction 
region (about 200 /im horizontal and 4 fim vertical) , is 
combined with the RMS of the transverse B flight length 
distribution to assign an uncertainty to the IP constraint. 

Simulation studies indicate that the vertexing proce- 
dure provides an unbiased estimate of zcp- The per- 
event estimate of the At error reflects the strong depen- 
dence of the Zcp resolution on the K'^ flight direction and 
the number of SVT layers traversed by its decay daugh- 
ters. For the 37% of events where both tracks include 
at least one hit in the inner three SVT layers (at radii 
from 3.2 cm to 5.4 cm), the mean At resolution is com- 
parable to that of decays for which the vertex is directly 
reconstructed from charged particles originating at the 
B decay point 13]. If both tracks have hits in the outer 
two SVT layers (at radii 9.1 cm to 14.4 cm) but one of the 
tracks has no hits in the inner three layers (^ 27% of the 
events) , the resolution is nearly two times worse. The re- 
maining events provide poor At measurements. For these 
events and for events with (TAt > 2.5 ps or ]At] > 20 ps, 
we do not include At information in the fit. However, 
we account for the contribution of these events in the 
measurement of C^o^o. 

We obtain the PDF for the time-dependence of signal 
decays from the convolution of Eq. ^ with a resolution 
function TZ(6t = At — Attruc, cta*)- The resolution func- 
tion is parameterized as the sum of a 'core' and a 'tail' 
Gaussian, each with a width and mean proportional to 
the reconstructed <TAt, and a third Gaussian centered 
at zero with a fixed width of 8ps 0]. We have veri- 
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fied in simulation that the parameters of TZ{dt, a^t) for 
B'^ K'^tt'^ decays are similar to those obtained from 
the -Bfiav sample, even though the distributions of a^t 
differ considerably. Therefore, we extract these parame- 
ters from a fit to the Bflav sample. We find that the At 
distribution of background candidates is well described 
by a delta function convolved with a resolution function 
with the same functional form as used for signal events. 
The parameters of the background function are deter- 
mined in the fit. 

To extract the CPV asymmetries we maximize the log- 
arithm of the likelihood function 

-(Ns+Ng) 

II [Nsfs4rs{x^,y^■, Sf, Cf) + iVfl/fle^Pfl(f„ y^■, q)] X 

igw/ At 

n [Nsil - fs)e'sr's{^,;Cf) + NBil - /b)6bPb(^«; a)] , 

iGw/o At 

where the second (third) factor on the right-hand 
side is the contribution from events with (without) 
At information. The probabilities Vs and Vb are 
products of PDFs for signal (S) and background 
(B) hypotheses evaluated for the measurements Xi = 
{toes, as, tag, tagging category} and yi = {At, aAt}- 
Along with the CPV asymmetries Sf and Cf, the fit ex- 
tracts the yields Ns and Nb , the fractions of events with 
At information fs and /s, and the parameters a which 
describe the background PDFs. 

Fitting the data sample of 4179 B" K*^n° can- 
didates, we find Ns = 122 ± 16 signal decays with 
Si^o^o = OAStoir ± 0.06 and C^^o^o = 0.40^|^;28 ± 0.09, 
where the uncertainties are statistical and systematic, 
respectively. The estimated number of signal decays is 
consistent with our measurement of the branching frac- 
tion The result for C^^o^o is consistent with a fit 
that does not employ At information. Fixing C^o^^o = 
we obtain 5*^0 ^^o — 0.41^0.48 ± 0.06. The evaluation of 
the systematic uncertainties is described below. 

Figure ^ shows the toes distributions for a signal- 
enhanced sample. The event selection is based on a like- 
lihood ratio R = VsHVb + Vs) calculated without the 
displayed observable. The dashed and solid curves in- 
dicate background and signal-plus-background contribu- 
tions, respectively, as obtained from the fit, but corrected 
for the selection on R. Figure |21 shows distributions of 
At for B^- and i3°-tagged events, and the asymmetry 
^^'0^0 (At) = [Nbo - Ngo] I \Nb° + ^Bo] as a function 
of At, also for a signal-enhanced sample. 

In order to investigate possible biases introduced in 
the CPV measurements by the IP-constrained vertex- 
ing technique, we examine i?*^ — > J/ijjK^ decays in data, 
where J/V' A*"*"/^" and J/V" e^e~. In these events 
we determine At in two ways: by fully reconstructing the 
B^ decay vertex using the trajectories of charged daugh- 
ters of the J/tp and the K'^ mesons, or by neglecting the 
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FIG. 1: Distribution of tties for events enhanced in signal de- 
cays. The dashed and solid curves represent the background 
and signal-plus-background contributions, respectively, as ob- 
tained from the maximum likelihood fit. 




FIG. 2: Distributions of At for events enhanced in signal de- 
cays with Stag tagged as (a) or (b) _B°, and (c) the asym- 
metry ^j^o^o(At). The dashed and solid curves represent the 
fitted background and signal-plus-background contributions, 
respectively, as obtained from the maximum likehhood fit. 
The asymmetry projection corresponds to approximately 36 
signal and 25 background events. 



J/V' contribution to the decay vertex and using the IP 
constraint and the K'^ trajectory only. This study shows 
that within statistical uncertainties the IP-constrained 
At measurement is unbiased with respect to the more 
established technique and that the obtained values of 
Sj^^i^o and Cj/^i^o are consistent. A similar study of 
B^ — > K'^n^ events, where the tt^ contribution to the 
decay vertex has been replaced by the IP constraint, 
yields Sko^± = 0.13 ± 0.19 and C^o^± = 0.06 ± 0.11, 



which is consistent with the expectation 5 



= 
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and our previous measurement of the charge asymme- 
try We also find that the B'^ hfetime measured in 
B° Kgir'^ decays and in IP-constrained -B" — *■ J/ipK'^ 
decays agrees with the world average [l^ . 

To quantify possible systematic effects we examine 
large samples of simulated B° K'^'^^ and _B" — > J/ip K'^ 
decays. We employ the difference in resolution function 
parameters extracted from these samples to evaluate un- 
certainties due to the use of the resolution function TZ 
extracted from the i?flav sample. We assign a systematic 
uncertainty of 0.03 on S'^^o^o and 0.02 on Cj^o^o due to 
the uncertainty in TZ. We compare fits to a large sample 
of simulated nominal and IP-constrained 5° J/ipK'^ 
events to account for any potential bias due to the ver- 
texing technique. This latter study yields the difference 
5Sj/^K° — 0.04, which we assign as the dominant sys- 
tematic uncertainty on 5*7^0 . We include a systematic 
uncertainty of 0.03 on S'/^o^o and 0.01 on C/^o^o to ac- 
count for a possible misalignment of the SVT. We con- 
sider large variations of the IP position and resolution, 
which we find to have negligible impact. We assign a 
systematic uncertainty of 0.09 to Cj^o^^o due to possible 
asymmetries in the rate of versus B'^ tags in back- 
ground events. Finally, we include a systematic uncer- 
tainty of 0.02 on both S'^j-o^o and Cj^o^o to account for 
imperfect knowledge of the PDFs used in the fit. 

In summary, we have performed a measurement of 
the time-dependent CPV asymmetries in K^ir'^ . 

These results supersede our previous measurement of 
C'xgu-o E3i which only relied on time-integrated observ- 
ables, and introduce the first measurement of S^o-kO- 
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